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ABSTRACT 

Computation  of  the  terms  of  the  vorticity  equation 
is  made  at  300  mbs  using  the  assumption  that  vertical  mo- 
tions are  negligible  and  approximating  relative  vorticity 
by  geostrophic  vorticity  computed  from  finite  differences. 
The  validity  of  the  computations  is  then  tested,  giving 
good  qualitative  agreement  but  poor  statistical  agree- 
ment,. The  isobaric  divergence  at  300  mbs  as  computed 
using  4£+tJ  =  Vp'V     is  compared  with  the  value  of 
divergence  as  obtained  from  the  vorticity  equation.  Again 
there  was  only  qualitative  agreement*  The  magnitude  of  the 
lateral  and  longitudinal  components  of  isobaric  divergence 
were  studied  and  it  was  found  that  the  lateral  component 
was  dominant  in  75   per  cent  of  the  cases  studied.  No  rela- 
tionship was  found  between  the  longitudinal  component  and  the 
total  divergence.  The  synoptic  aspects  of  the  divergence 
fields  at  300  mbs  are  discussed  with  respect  to  the  500-mb 
surface. 

The  writer  wishes  to  express  his  appreciation  for  the 
guidance  and  encouragement  given  him  by  Professor  F.  L.  Martin 
of  the  U.S.  Naval  Postgraduate  School  in  this  investigation. 
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1.  Introduction. 

Several  studies  of  divergence  in  the  atmosphere  have 
been  carried  out  by  different  methods  and  by  different  inves- 
tigators.  Some  of  the  most  notable  of  these  are  Fleagle  \2~]   , 
Arnason  [l]  ,  and  Landers  [5]  .  Their  investigations  have 
involved  a  three  dimensional  study  of  the  distribution  of 
divergence  in  the  atmosphere.  Landers  [5]  computed  horizon- 
tal divergence  using  constant  level  charts  and  a  "del  com- 
puter" which  consists  of  computing  divergence  from  the  wind 
vectors  at  the  vertices  of  an  equilateral  triangle  with  an 
altitude  of  2&°   of  latitude.  Arnason  [l]  computed  the  inte- 
grated divergence  over  a  12  hour  period  in  isobaric  layers 
with  the  use  of  vertical  velocities,  the  local  rate  of 
change  of  thickness  and  the  temperature  advection.   Fleagle 
computed  horizontal  velocity  divergence  and  mass  divergence 
up  to  about  150  mbs  and  found  a  level  of  non-divergence  at 
about  600  mbs. 

This  study  is  carried  out  at  300  mbs  only  and  no  attempt 
is  made  to  extend  it  to  three  dimensions.   The  main  objective 
of  this  study  was  to  gain  some  practical  knowledge  of  the 
behavior  of  divergence  patterns  and  to  apply  the  results  to 
the  prognoses  of  the  500-mb  chart.   Computations  of  the  di- 
vergence and  the  individual  rate  of  change  of  the  absolute 
vorticity  with  time  are  made,  and  their  relationship  tested. 
Also  a  study  of  the  relationship  between  the  longitudinal 
component  of  the  divergence  and  the  isobaric  divergence  is 
made. 


2.  Source   of  Data. 

A  series  of  300-mb  maps,  carefully  analyzed  by  the 
staff  of  the  U.S.  Naval  Postgraduate  School,  was  used  in 
making  the  computations.  The  analyses  employed  were  those 
of  1500Z,  15  November,  0300Z  and  1500Z,  16  November  1955. 
On  the  first  chart  there  is  a  trough  and  a  closed  low 
at  about  112°  W  and  a  ridge  at  about  8*+°  W  with  consider- 
able meridional  flow.  At  0300Z  on  the  16th  the  trough  and 
ridge  have  moved  about  5°  eastward,  retaining  about  the 
same  intensity..  By  1500Z  on  the  16th. the  amplitude  of  the 
trough-ridge  system  has  decreased  considerably  and  moved 
to  98°  W  and  72°  W  longitude  so  that  in  between  the  trough 
and  ridge  the  flow  is  becoming  zonal..  Computations,  which 
will  be  described  in  the  following  sections,  were  carried 
out  between  this  trough  and  ridge.  In  addition  some  com- 
putations were  made  behind  the  trough  on  the  0300Z,  16 
November  chart..  A  total  of  138  computations  of  all  five 
terms  of  Equation  (7)  below  were  made. 

3.  Mathematical  Analysis. 

The  vorticity  equation  in  isobaric  coordinates  is: 

CD    4^  =  -mvw-  u- vPu)  xi^H 

dt  "  ^  p 

This  is  not  in  the  form  in  which  computations  from  a 
synoptic  chart  can  be  readily  made.   Since  vertical  veloc- 
ities are  normally  small,  the  second  term  on  the  right  side 


of  the  equation  may  be  ignored..  Equation  (1)  then  becomes 
the  approximation  to  the  vorticity  equation  as  it  is  common- 
ly made  in  meteorology. 

Making  the  geostrophic  approximation,  T=  3.  VLz   »  and 
noting  that  the  coriolis  parameter  does  not  change  with 
time,  Equation  (2)  becomes: 

In  order  to  compute  the  individual  rate  of  change  of 
vorticity  it  would  be  necessary  to  construct  isobaric  tro- 
jectories  of  the  air  parcels.  However,  such  computations 
may  be  avoided  by  making  use  of  the  fact  that: 

Dt    }t    ~H  ~P 

Equation  (3)  then  becomes" 

J 

Since  £  is  quasi-constant  compared  to  the  horizontal  var- 

f 
iation  of  \7  z  ?  it  may  be  brought  outside  of  the  differential 

operator  in  the  second  term  of  the  left  side  of  Equation  (5). 

In  order  to  get  v7  z   in  a  form  in  which  computations  can  be 

made,  finite  differences  may  be  used  to  give: 

(6)   V\   *  ±p- 


Combining  Equations  (5)  and  (6)  and  multiplying  by  fd2  gives: 

%  g 

(7)  ^-z)  ♦  &•&<*-*)   +  y?"  =  -[(*-*) +^](7p-V) 


The  equation  is  now  in  a  form  in  which  computations  were 

made.  The  left  hand  side  of  this  equation  now  represents 

fd  D  J«.   and  the  right  hand  side  represents;  -fd  ^(Vp-V/J. 

g  ET  g 

By  measuring  wind  velocities  in  feet  per  second  and  grid 

distances-  in  feet,  the  units  of  this  equation  can  be  shown 

to  be  feet  per  second. 


hi   Computational  Procedure. 

The  quantity  z  may  be  readily  computed  using  a  grid 
as  shown  in  Figure  1. 


Figure  1 

Grid  for  computing  space  mean. 

<> 

The  value  of  2d  was  chosen  to  be  5°  of  latitude.  This 
value  ±m  about  half  of  the  length  of  the  short-wave  trough 
so  that  in  constructing  a  "Z -chart,  the  short-wave  troughs 
tend  to  be  removed,  but  the  long-wave  troughs  are  unaf- 
fected. 
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The  construction  of  the  z-chart  where  z  is'  z^-t-  Z2+  Z3+  zi+ 

Tf 

(Figure  1),  was  accomplished  by  displacing  two  z-charts  10 
degrees  of  latitude  laterally  and  adding  graphically  giving 
z^+  Z3.  The  two  maps  were  then  displaced  10  degrees  of 
latitude  northward  and  added  giving  z24Zi+.  The  Z1+Z3 
contours  are  then  added  to  the  z2+zif  contours  giving  Vz. 
By  dividing  by  h   we  then  have  a  space  mean  map  from  which 
the  actual  map  may  be  graphically  subtracted  giving  iso- 
lines-  of  Oz-z).  The  local  rate  of  change  of  (z-z)  was 
approximated  by  the  average  rate  of  change  in  the  quantity 
over  a  12-hour  interval,  as  measured  at  any  given  station. 
The  advection  term  [term  4b )J  was  computed  by  measur- 
ing (1)  the  gradient  of  (z-z)  across  selected  stations 
over  a  distance  of  5°  of  latitude,  (2)  the  wind  velocity 
in  feet  per  second,  and  (3)  the  angle  p    between  the  wind 
vector  and  the  gradient  vector..  The  magnitude  of  term  (b) 
is  then  |Vl  |VP(Z-z)|  \c°s  ¥\  .     A  diagram  of  the  vectors 
involved  is  shown  in  Figure  2. 


Vf(*-z) 


Figure  2 
> 

Diagram  of  Angle  ip 

All  factors  of  term  (c)  with  the  exception  of  v  are 
constant  for  a  given  latitude  and  can  be  computed  for  each 
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latitude  circle.  Taking  the  appropriate  value  of  fd^  ft 

g 
and  multiplying  "by   the  northward  component  of  the  wind 

velocity  gives  the  final  term  in  the  vorticity  computation. 
Now  by  equation  (7)  the  sum  of  terms  (a)-,  (b)  and  (c)  rep- 
resent the  quantity  fd2  D  &  subject  to  the  accuracy  of  the 

g  Dt 
approximations  made  in  neglecting  vertical  motions  and  in 

computing  _^_  ("z-z) . 
it 
Several  methods  have  been  devised  for  computing  iso- 

baric  divergence1  |3>  6  J  .   It  was  decided  to  use  the  fol- 
lowing equation  in  the  computation. 

A  square  grid  2%   degrees  of  latitude  on  a  side  was  used. 

An  average  u  was  computed  on  each  side  of  the  square  and 

an  average  v  at  the  top  and  bottom.   As  shown  in  Figure  3> 

the  isobaric  divergence  is  then  u2-U]_  v2-v^ 

d      (I 


Hi 


V, 


4i 


d»2| 


i  • 


I* J'*i' >j 

Figure  3 
Grid  for  computing  divergence. 


Since  the  grid  is  square  no  spherical  correction  is 
necessary,  and  grids  can  be  oriented  along  the  major  flow. 
If  the  values  of  isobaric  divergence  are  then  multiplied 


by  tTwm  (d)  of  Equation  (7)  the  value,  should  be  the  same 

as  the  sum  of  terras  (a),  (b)  and  (c).  The  divergence  was 

also  computed  by  dividing  the  sum  of  terms  (a),  (b)  and  (c) 

by  term  (d)  of  Equation  (7).   For  the  discussion  in  this 

paper  this  value  of  isobaric  divergence  will  be  referred 

to  as  effective  divergence. 

Isolines  of  the  computed  values  of  fd2  D  fo.  are 

g  Dt 

plotted  in  Figures  h,   8  and  12.  The  isolines  of 

-fd2  jriT Jyr-v)     are  drawn  in  Figures  5,  9  and  13.  The  field 

g 
of  effective  divergence  is  plotted  in  Figures  6,  10  and  Ih, 

and  isobaric  divergence  is  plotted  in  Figures  7?  11  and  15 '. 

The  labeled  values  of  isolines  of  fd  D  ffy   and 

2      x  g  Dt 

-.£1  f«(¥f"Y/   are  times  10  J  and  the  labeled  values  of  the 

g  -5" 

fields  of  divergence  are  times  10  .   In  plotting  the  values 

of  isobaric  divergence  in  Figure  7  the  values  in  four  adja- 
cent squares  were  averaged  and  then  plotted  in  the  common 
point  of  the  four  squares..  This  was  done  in  an  effort  to 
smooth  the  divergence  field  but  was  not  done  in  Figures  11 
and  15  because  the  averaging  could  mask  the  effects  of  short 
waves. 

5.  The  Results  of  the  Computations. 

The  individual  terms  on  the  right  hand  side  of  Equation 

-3 

7  have  order  of  magnitude    10  '  feet  per  second.  A  typi- 
cal value  of  the  change  of  Cz-z)  over  a  twelve  hour  period 
is -100  feet.  When  divided  by  12  hours  or  V.32  x  10  seconds, 


-^x  (z-z)  has  a  value  of  about  t  2.3  x  K^^feet  per  second. 

The  extreme  computed  value  of  -i—  ("2-z)  was  11.6  x  10   feet 
per  second.  Measuring  \V9(-t;-z)\  over  a  distance  of  5  of  lati- 
tude (2^-°  each  side  of  the  station)  gives  values  of  the  order 

-5 
of  5.0  x  10  .  With  a  wind  of  120  feet  per  second  and  an 

* 

angle  <p   whose  cosine  is  0.3 »  a  typical  value  of  term  (b) 

is  1.8  x  10   feet  per  second.  The  maximum  magnitude  of 

_3 
this  term  was  20.5  x  10   feet  per  second.   The  quantity 

£L  (3    ranges  from  11.1+  x  i0"6  at  30°  and  60°  N  to 
13.18  x  10"  '  at  *+5°  N  where  it  is  maximum.  When  multiplied 
by  a  typical  value  of  100  feet  per  second  for  v,  one  obtains 
a  value  for  term  (c)  of  1.1  x  10"^  to  1.3  x  10~^  feet  per 
second.   The  maximum  value  of  term  (c)  was  I-87  x  10"^  feet 
per  second.  Although  term  (c)  is  of  the  same  order  of  magni- 
tude as  terms  (a)  and  (b) ,  it  is  somewhat  smaller  in  general 
and  terms  (a)  and  (b)  give  the  largest  contribution  to  the 
sum  of  the  three  terms  of  the  left  side  of  Equation  (7) . 
The  value  of  f2d2  is  constant  for  a  given  latitude 

S.  O  r\  O 

circle  and  ranges  from  1.37  x  10  feet  at  30,  N  to  3*66  x  10 
feet  at  55°   N,  while  (z-z)  ranges  from  -300  feet  to  300  feet. 
Term  (d)  was  positive  in  all  the  computations  except  one,,  in 
which  it  was  a  small  negative  value.   Of  course  all  compu- 
tations were  carried  out  in  the  broad  scale  flow  between 
trough  and  ridge.  Term  (d)  could  very  well  be  a  negative 
value  in  a  closed  high  pressure  system  where  values  of 
(z-z)  might  be  as  great  as  -UO0  feet. 

Values  of  isobaric  divergence  were  in  general  smaller 
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than  those  of  effective  divergence.  However  they  agree 
fairly  well  qualitatively-  The  magnitude  of  the  effective 
divergence  agrees  very  well  with  that  computed  by  Landers  [_5]* 
The  average  absolute  value  of  the  effective  divergence  on 
the  three  charts  was  3.3*+,  3.70  and  1.22  all  times  10" ^  per 
second.   A  plot  of  divergence  values  versus  height  was  made 
by  Landers  Qj?]  and  shows  an  increase  of  0.6  x  10"^  per 
second  in  the  mean  absolute  value  of  horizontal  divergence 
from  20,000  to  25,000  feet.  Extrapolation  of  his  curve  to 
30,000  feet  would  give  a  value  of  2.7  x  10   per  second,  and 
the  mean  absolute  value  of  effective  divergence  of  all  the 
computation  made  in  this  study  is  2.78  x  10""?  per  second. 

The  values  of  isobaric  divergence  have  a  maximum  magnitude 

-5 
of  about  3.5  x  10   per  second  while  the  mean  absolute 

value  is  about  1.0  x  10   per  second  which  agrees  better 

with  the  values  of  divergence  as  computed  by  Arnason  £lj  . 

In  the  following  section  a  statistical  approach  is  made  in 

comparing  the  two  fields  of  divergence  and  the  fields  of 

fd£  DJL  and  -fd2  -Mtyv) 
g   Dt        g 

6..  Statistical  Studies. 

Of  the  patterns  of  fd2  D  &  and  -fd2  ^p  .y) 
Figures  12  and  13  seem  to  give  the  best  agreement.   To  test 
this  agreement  a  correlation  coefficient  of  the  two  quantities 
was  found  using  values  at  35  points  in  each  figure.   This 
correlation  coefficient  r  is  given  by  the  following  equation 
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from  Hoel  (>]. 

(9)  r  -  cov  (xy) 

sx  sy 

Cov(xy)  is  the  covariance  of  x  and  y,  the  stochastic  variables, 

and  Sx  and  Sy  are  the  respective  standard  deviations.   Cov(xy) 

and  the  standard  deviations  were  computed  from  the  following 
equations  gi^en  in  [V]. 

(10)  cov(xfc)  -  -^  [klty  -  ^  *  jj 

to)  s§  =7j[**E-*fe<] 

(12)  Sf  ,  x[<^-^f] 

The  mean  value  of  fd2  D  &.  was  -2.06  x  10^  feet  per 

S  Dt  \  -3 

second  and  the  mean  value  of  -fd£  1*,{V?'V  )  was  0.1  x  10 

6  -3 

feet  per  second.  The  standard  deviations  were  6,38  x  10 

feet  per  second  and  the  3.52  x  10"  3  feet  per  second  respect- 
ively. The  correlation  coefficient  was  found  to  be  0.31. 
With  90  per  cent  confidence  limits -the  critical  value  of 
r  is  0.29  so  the  correlation  coefficient  of  0.31  is  signif- 
icant at  this  level  of  confidence.  However,  the  value  of 
r  is  not  large  enough  to  deduce  more  than  that  the  two 
quantities  are  related.   Some  of  the  reasons  for  this  lack 
of  agreement  will  be  discussed  later. 

The  same  procedure  as  described  in  the  previous  para- 
graphs was  used  to  compute  a  correlation  coefficient  between 
Figures  10  and  11  which  are  patterns  of  effective  divergence 

10 


and  isobaric  divergence.   It  vas  believed  that  perhaps  better 
agreement  could  be  achieved  between  these  two  patterns  since 

in  the  process  of  dividing  D  j'q  by -f  some  of  the  errors  in 

Dt 
using  the  geostrophic  approximation  might  be  minimized.  How- 
ever, the  correlation  coefficient  between  Figures  10  and  11 
was  only  0.1^  which  is  significant  at  68  per  cent  confidence 
limits  with  a  sample  size  of  52. 

The  approximation  of  _£_  (£-z)  by  a  12-hour  average 
would  tend  to  cause  a  computed  effective  divergence  center 
to  be  east  of  its  instantaneous  position  if  the  future  12- 
hour  Cz-z)  map  is  used  and  west  of  its  instantaneous  posi- 
tion if  the  previous  12-hour  ("z*-z)  map  is  used.   The  plotted 
values  of  effective  divergence  and  isobaric  divergence  were 
examined  in  an  attempt  to  find  some  consistent  displacement 
which  would  give  a  higher  correlation  coefficient.   It  was 
believed  that  the  correlation  coefficient  could  be  improved 
by  shifting  the  divergence  center  over  Indiana  in  Figure  10 
into  coincidence  with  the  divergence  center  over  Michigan 
in  Figure  11  while  the  convergence  centers  over  Iowa  remained 
in  coincidence.  Actually,  however,  there  was  no  improvement 
in  the  correlation  coefficient.  A  value  of  0.07  was  com- 
puted for  r  with  a  sample  size  of  31  which  is  not  significant 
at  any  reasonable  level  of  confidence. 

Sutclif  fe  [?~\   has  stated  that  the  vortex- tube  term  in  • 
the  vorticity  equation  tends  to  be  offset  by  -  ¥(¥p'  Y/ 
so  that  a  better  approximation  to  the  vorticity  equation 
would  be : 
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(13)      JD2*  -  -f(?  -v/  ) 
Dt      \-r  -  / 


A  correlation  coefficient  was  computed  using  this 
approximation  and  the  computations  from  the  1500Z,  16 
November  chart.  Using  Equations  (9),  (10),  (11)  and  (12) 
the  correlation  coefficient  was  computed  to  be  -0.201!-  which 
is  significant  at  75  per  cent  confidence  limits.   The  fact 
that  the  correlation  coefficient  is  negative  seems  to  indi- 
cate that  ignoring  j\?r"Y)must  overcompensate  for  the  vortex- 
tube  term. 

Using  the  equation  for  isobaric  divergence  in  natural 
coordinates: 

it  is  seen  that  there  are  two  components  of  the  total  di- 
vergence, a  longitudinal  component  along  a  streamline  and 
a  lateral  component  normal  to  the  streamline.   It  was 
also  of  interest  to  determine  which  of  these  components 
gives  the  greatest  contribution  to  the  total  isobaric  di- 
vergence. The  rate  of  change  of  wind  velocity  along  a 
streamline  is  most  easily  measured  and  was  found  at  all  the 
stations  at  which  vorticity  computations  were  made.  This 
was  done  by  measuring  the  change  in  wind  speed  along  a 
streamline  for  a  distance  5°  of  latitude  across  each  sta- 
tion.  If  these  values  are  compared  with  the  effective 
divergence,  the  difference  will  be  the  effective  lateral 

component.  In  most  cases  the  lateral  and  longitudinal 
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terms  are  of  opposite  sign.   In  103  out  of  138  computa- 
tions (approximately  7%)    the  lateral  term  was  the  domin- 
ant term.   Therefore,  most  of  the  time  the  sign  of  the 
effective  divergence  is  determined  by  the  effective 
lateral  component. 

It  would  be  helpful  if  some  relationship  could  be 
derived  between  effective  divergence  and  the  longitudinal 
component  since  the  longitudinal  component  can  be  easily 
measured.  'The  total  divergence  could  then  be  determined 
simply  by  measuring  ^L    .  Using  Equations  (9),  (10),  (11) 

dS 

and  (12)  to  determine  a  correlation  coefficient  between 
-^-  and  the  effective  divergence  gave  a  coefficient  of 
O.O87  which  is  not  significant  at  any  reasonable  level 
of  confidence.   Therefore,  very  little  about  the  magnitude 
of  the  effective  divergence  can  be  deduced  from  the  magni- 
tude of  the  longitudinal  divergence. 

7.  Sources  of  Error. 

The  discrepancy  which  exists  between  computed  values  of 

2  2 

fcL  PJL  and  -fd_  \(a,\y)      and  between  effective  and  iso- 

baric  divergence  is  due  partially  to  the  fact  that  the  vortex- 
tube  term  and  the  vertical  advection  of  vorticity  have  been 
ignored.  However,  the  discrepancy  is  greater  than  the  values 
to  be  expected  from  the  terms  that  have  been  ignored.  Due  to 
the  lack  of  actual  wind  reports  at  300  mbs,  it  was  necessary 
to  interpolate  between  the  actual  wind  observations.  Error 
could  have  been  introduced  here*  Also  a  10°  error  in  the 

13 


observed  wind  direction  and  a  5-knot  error  in  wind  speed 
may  cause  a  maximum  error  of  3.0  x  10"  ?  per  second  in  the 
computed  isobaric  divergence  while  the  same  error  in  wind 
direction  and  speed  will  cause  only  an  error  of  0.3*+  x  10"^ 
in  the  effective  divergence.  However,  there  is  another 
approximation  in  the  use  of  the  vorticity  equation.  This 
was  the  determination  of  ~-  (z-z)  as  a  12-hour  average 
using  two  consecutive  (z-z) -charts.  The  magnitude  of  the 
errors  due  to  this  approximation  is  unknown. 

8.  Synoptic  Aspects. 

It  was  difficult  to  determine  the  track  of  the  diver- 
gence centers  as  many  of  them  disappeared  or  were  out  of  the 
•area  of  computations  in  the  12-hour  intervals*  For  example, 
in  the  12-hour  interval  from  Figure  6  to  Figure  10,  three 
divergence  centers  have  appeared  where  there  was  previously 
only  one. 

The  effective  divergence  patterns  give  good  qualitative 
agreement  with  the  synoptic  situation..  Fleagle  [2]  has  found 
the  500-mb  level  to  be  approximately  a  level  of  non-divergence 
and  the  300-mb  surface  is  near  the  level  of  maximum  divergence. 
Landers  [5J  also  found  the  300-mb  level  to  be  a  level  of  maxi- 
mum divergence  and  the  500-mb  level  to  be  a  level  of  least 
divergence.  Arnason's  results  £l]  did  not  show  such  a  distri- 
bution of  divergence,  but  it  is  believed  that  the  isobaric 
divergence  as  computed  at  300  mbs  may  give  a  good  indication 
of  rises  and  falls  of  the  500-mb  surface.   In  Figure  6  the 


areas  of  effective  divergence  between  trough  and  ridge  are 
relatively  large  and  clear  cut.  With  this  divergence 
pattern  aloft  we  would  expect  the  trough  to  continue  to 
move  eastward,  which  actually  occurs.  Figure  7  does  not 
give  this  same  agreement  as  there  is  a  convergence  center 
centered  over  station  53^  which  is  behind  the  ridge.  This 
would  indicate  retrogression  of  the  ridge  which  does  not 
occur.   It  may  also  be  noted  in  Figure  6  that  as  the  strong 
winds  east  of  the  ridge  at  8^°  W  enter  the  cyclonically 
curved  isohypses  of  weaker  gradient  there  is  a  change  from 
convergence  to  divergence.   On  the  left  hand  side  of  Figure 
6  there  appears  to  be  a  perturbation  leaving  the  long-wave 
trough  as  indicated  by  the  divergence  center  over  Colorado. 
As  yet  there  is  no  actual  short-wave  trough  indicated  at 
500  mbs.   On  the  next  500-mb  chart,  0300Z  16  November,  a 
short-wave  trough  has  appeared  east  of  the  long- wave  trough. 
Figures  10  and  11  indicate  closed  divergence  centers  ahead 
of  the  short-wave  trough  and  closed  convergence  centers 
behind  it.  With  the  short-wave  trough  in  the  broad  scale 
flow,,  the  pattern  of  divergence  is  much  more  complicated  and 
the  gradients  are  much  greater.   On  the  West  Coast  in  Figures 
10  and  11  there  is  a  center  of  convergence  located  ahead  of 
a  short-wave  trough.   Twelve  hours  later  the  trough  is  prac- 
tically non-existant.  The  divergence  center  located  over 
Oregon  may  be  associated  with  the  trough  that  appears  off 
the  West  Coast  on  the  1500Z  16  November  chart.   In  Figures 
Ik   and  15  the  short-wave  trough  of  Figures  10  and  11  fcas 
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moved  into  the  ridge  and  weakened.   Consequently  the  divergence 

and  convergence  centers  that  were  associated  with  it  are 

not  ,as  clear  cut.   The  patterns  have  become  elongated 

along  the  contours  with  considerable  convergence  east  of  the 

trough.   Therefore,  the  trough  may  be  expected  to  decrease  in 

amplitude  and  the  flow  to  become  more  zonal,  which  is  verified 

on  the  0300Z,  17  November  500-mb  chart* 

Due  to  the  fact  t^iat  the  surface  maps  are  two  and  a 
half  hours  earlier  than  the  upper-air  charts,  it  is  diffi- 
cult to  compare  divergence  patterns  with  the  surface  maps. 
However,,  in  the  few  cases  studied  the  low-pressure  centers 
at  the  surface  move  toward  the  centers  of  divergence  at 
300  mbs*  It  is  also  interesting  to  note  that  areas  of  di- 
vergence tend  to  overlay  the  fronts  at  the  surface.. 

9.  Summary  and  Conclusions. 

This  investigation  of  the.  isobaric  divergence  patterns 
at  300  mbs  has  shown  that  the  technique  of  computing  diver- 
gence from  the  rate  of  change  of  absolute  vorticity  gives 
very  good  agreement  in  magnitude  with  the  computations  of 
isobaric  divergence  made  by  Landers  [5}   .  The  computations 
made  in  this  investigation  indicate  that  a  mean  absolute 
value  of  divergence  at  300  mbs  is  about  2.8  x  10"*?  per  second* 
This  value  of  effective  divergence  acting  over  a  column  of 
air  12,000  feet  thick  for  a  period  of  12  hours  will  cause  a 
fall  in  the  500-mb  surface  of  about  200  feet.   The  values  of 
effective  divergence  plotted  in  Figures  6,  10  and  l*f  give 

good  indications  of  future  areas  of  rise  and  fall  at  the 
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500-mb  surface. 

The  computations  of  isobaric  divergence  using  Equation 
(8)  gave  a  smaller  mean  value,  which  gives  good  agreement 
with  the  magnitudes  computed  by  Arnason  |_lj  .  Except  in 
Figure  7 »  these  values  of  divergence  also  give  good  indi- 
cations of  future  developments  on  the  ^OO-mb  surface. 

There  were  rather  poor  statistical  correlation  co- 
efficients between  computed  values -of  fd2  D^a   and 

2  S  W 

-£L-  j^ (yp  •  V  ),   and  also  between  effective  and  isobaric  diver- 

g       "" 
gence..  This  poor  agreement  is  believed  due  largely  to  errors 

in,  and  insufficient,  wind  reports  at  300  mbs»  Also  the 

12-hour  approximation  of  -ir(z-z)  introduces  error. 

The  lateral  component  of  the  effective  divergence  was 

the  dominant  component  in  75%   of  the  computations  made* 

There  was  found  to  be  no  relationship  between  the  effective 

divergence  and  either  the  lateral  or  longitudinal  component 

of  divergence.. 
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